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Abstract 

We investigate a recent proposal to construct chiral gauge theories on the lattice using 
domain wall fermions. We restrict ourselves to the finite volume case, in which two 
domain walls are present, with modes of opposite chirality on each of them. We couple 
the chiral fermions on only one of the domain walls to a gauge field. In order to preserve 
gauge invariance, we have to add a scalar field, which gives rise to additional light 
mirror fermion and scalar modes. We argue that in an anomaly free model these extra 
modes would decouple if our model possesses a so-called strong coupling symmetric 
phase. However, our numerical results indicate that such a phase most probably does 
not exist. 
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1 Introduction 



The lattice provides a first-principles regularization of quantum field theories, which 
allows us to explore the nonperturbative properties of a model and for vectorlike theories, 
such as QCD, it has proven to be very successful. Since the full Standard Model is a 
chiral gauge theory, it is natural to attempt a construction of chiral gauge models on 
the lattice as well. 

As is well known, however, on the lattice one is confronted with "species doubling" 
m, 0, i.e. the phenomenon that a single Weyl fermion field on the lattice leads to an 
equal number of left and right handed fermions in the continuum limit. When coupled 
to a gauge field, all doublers transform in the same representation of the gauge group 
which prevents an easy construction of chiral gauge theories. For nonchiral models there 
are two well tested ways of dealing with the species doublers: they can be decoupled 
with a momentum dependent mass term as in Wilson's method, or they can be used as 
physical degrees of freedom as in the staggered fermion method. However, since these 
methods violate chiral symmetry, a straightforward extension to chiral gauge theories 
would clash with gauge invariance. Proposals for chiral gauge theories on the lattice 
include generalizations of Wilson's method [3-6] and of the staggered method 0. There 
are proposals that try to avoid coupling the doublers [Q, and approaches that start 
from a gauge- fixed continuum action ||^, 0. There are also proposals for more radical 
departures from the usual lattice fermion prescriptions [11-13]. For a recent review, see 
ref. [0. 

The domain wall fermion approach suggested in ref. [|12] falls into the last group 
and has attracted a lot of attention recently [15-22] . In the domain wall model an extra 
dimension is added to our four dimensional world. In this five dimensional world the 
model is vectorlike and the fermion doublers can be removed using Wilson's method 
without breaking gauge invariance. The reduction to a four dimensional world with a 
chiral fermion is made by giving the fermions a mass term which flips sign across a four 
dimensional domain wall. It has been shown [0 that the lattice Wilson-Dirac operator 
with such a mass term has a chiral zeromode, which is bound to the domain wall. This 
fermion remains massless and localized at the domain wall for (four-)momenta below a 



critical momentum [|l2i |T5|, [T^]. On a finite lattice the (periodic) boundary conditions 
lead to a second anti-domain wall with a chiral fermion of opposite handedness. 

Every lattice model for a chiral gauge theory has to produce the appropriate anomaly 
structure of the target continuum theory. The domain wall model has the potential 
to solve this problem elegantly with the help of the extra dimension. The starting 
five dimensional model is vectorlike and hence the gauge current is anomaly free, 
X]^=i dfj^J^ = 0. However, the four dimensional current restricted to the domain wall is 
clearly not conserved Y^t=id^J^ = —d^J^, and its divergence reproduces the expected 



anomaly when computed for weak external gauge fields [12, 15, M. J5 takes the form 



of a Goldstone-Wilczek current with a nonzero derivative across the domain wall, as 
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was demonstrated some time ago in the continuum in ref. [Q. On a finite lattice this 



Goldstone-Wilczek current transports charge from the domain wall to the anti-domain 
wall, ensuring charge conservation in the five dimensional theory. The same mechanism 
should also yield the correct four dimensional global anomaly structure. 

In the work referred to above the domain wall fermions are coupled to fixed smooth 
external gauge fields. Here it is not important that the chiral fermions at both the 
domain and anti-domain walls couple to the gauge field, because we can single out 
one of the domain walls by hand. With dynamical gauge fields, however, the crucial 
requirement is that only a single domain wall fermion couples to the gauge field. If this 
can be achieved the second domain wall can be ignored and we are left with an interacting 
chiral fermion in a four dimensional world located at the domain wall, assuming that 
the fermions are in an anomaly free representation. 

In the original proposal it was hoped that the communication between the two do- 
main walls could be prevented by modifying the gauge interactions in the fifth dimension. 
However, it seems likely that this approach does not lead to the desired decoupling (see 
also ref. |2^) and here we follow instead the suggestion made in ref. In this ap- 

proach the gauge fields are coupled only in a restricted region around one of the domain 
walls, where the size of this region, which we will call the waveguide, should be at least 
as large as the support of the wave function of the domain wall zeromode. However, as 
will be discussed in much more detail below, the requirement of gauge invariance leads 
to the introduction of an extra scalar field at the boundaries of the waveguide. This 
scalar field screens the gauge charge of the fermions at the waveguide boundary and 
allows for interactions between these charged fermions and the neutral ones outside the 
waveguide. This leads to Yukawa couplings located at the waveguide boundary which 
give rise to additional light fermion modes at the waveguide boundary. This is most 
easily seen for zero Yukawa coupling, because then the waveguide region decouples from 
the anti-domain wall region and five dimensional charge conservation is now ensured by 
the new zeromodes at the boundary. 

The aim is then to decouple the fermion at the waveguide boundary and maintain at 
the same time the chiral zeromode at the domain wall. Because we have introduced a 
scalar field coupled to the fermions through a Yukawa interaction, we may hope for a rich 
phase structure of the model, similar to that found in other two and four dimensional 
Yukawa models on the lattice. In particular, one expects that one can drive the system 
into a symmetric phase, with vanishing scalar field vacuum expectation value f , and a 
spontaneously broken phase, with v > 0. For small values of the Yukawa coupling one 
then expects the fermions at the waveguide to follow the perturbative relation mp ocv 
(with f = in the symmetric phase). This means that these fermions remain light and 
appear in the low energy spectrum. However, at large values of the Yukawa coupling, the 
interaction of the fermion and scalar fields might become so strong that only a bound 
state fermion exists with a mass of the order of the cutoff. Such a strong coupling 
behavior has been established in various Yukawa models on the lattice [25-27]. If the 
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four dimensional model is anomaly free and if we choose the waveguide boundary with 
the scalar field far enough from the domain wall, we could hope to take the waveguide 
fermions into a strong coupling symmetric phase, without affecting the chiral mode at 
the domain wall. Then the fermions at the waveguide boundary would decouple from 
the low energy physics, leaving only the chiral zeromodes at the domain wall coupled to 
the gauge field. 

The crucial question we will investigate in this paper is therefore whether such a 
strong symmetric phase exists in the domain wall model. We will present evidence 
based on analytical considerations and numerical results, which leads us to conclude 
that an appropriate strong coupling phase most probably does not occur. 

The paper is organized as follows. In sect. 2 we review free domain wall fermions, 
discuss the coupling to the gauge field and the need to introduce the extra scalar field. 
We close the section with a sketch of the phase diagram we would hope to find for our 
model. In sect. 3 we rewrite the fermion action in a mirror- fermion form, such that we 
can distinguish the light modes from the heavy ones. In the next section we present 
results for fermion masses, concentrating on the results for the boundary fermion. In 
sect. 5 we continue our search for a strong coupling phase using the eigenvalue spectra 
of the fermion matrix in a simplified model, in which all heavy modes are discarded. 
Sect. 6 contains a brief discussion of alternative ways to couple the gauge field and in 
sect. 7 we present our conclusions. 



2 Domain wall fermions coupled to gauge fields 
2.1 Resume of free domain wall fermions 

Let us start our discussion with a short resume of free domain wall fermions. Consider an 
odd dimensional lattice of size L'^Lg, with d = 2n, Lg the extent in the extra dimension 
and lattice sites labeled by (x, s), (x = (xi, ■ ■ ■ , x^)). The action for free domain wall 
fermions [0 can be written as, 

s \xy 

- hT.l^lir-j,)%-''+X:'\r + ^,)^:-2rr^n]), (2.1) 



where ^ and w are the Dirac operator and Wilson term with Wilson parameter r on 
the even d dimensional lattice, 

d 
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d 

The 5x^y is the Kronecker delta and we use lattice units a = 1. We shall choose the 
Wilson parameter r = 1 for convenience and for the domain wall mass, denoted by m^, 
we choose a periodic step function of the form (with Ls even), 

= -mo s = 2,- ■ ■ ,Ls/2, 
m' = s = l,L,/2 + l, (2.3) 
m'* = +mo s = Ls/2 + 2,- ■ ■ ,Ls. 

With periodic boundary conditions in s, the emergence of an anti-domain wall is in- 
evitable. It will often be convenient to think of s as a flavor label, rather than an extra 
space-time coordinate. 

This model posesses two chiral zeromodes with the property that the mode bound 
to the domain wall at s = 1 is left handed (75 = — 1) and the mode bound to the other 
domain wall at s = Ls/2 + 1 is right handed (75 = +1) ||12[. The wave functions for both 
modes have the form of plane waves in the 2n- dimensional space and decay exponentially 
in s, away from the domain walls. These chiral zeromodes exist for plane wave momenta 
below some critical momentum which depends on the ratio mo/r. For different values 
of mo/r the zeromode spectrum can change substantially. For < mo/r < 2 one has only 
one chiral zeromode at each domain wall. For increasing values of mg/r this zeromode 
becomes less localized and disappears at mg/r = 2. At this point new zeromodes with 
opposite chirality are provided by the species doublers which are located at different 
corners of the Brillouin zone [|I7|, |TB|. Throughout the paper we will take mo/r 1 
and hence we will have only one chiral zeromode at the domain wall with exponentially 
small overlap with the zeromode at the anti-domain wall. In this case, chiral modes 
exist for momenta p below a critical momentum, \p\ < p^, with p^ = 2 1]^(l ~" cos {p^)) 
and p1 = A — 2mo/r. Note that p^ ^ p^ for small momenta. 

2.2 Coupling to gauge fields 

Since the left and right handed zeromode components of the fermion field now are 
separated in s space, one can attempt to couple these two components in different ways 
to a gauge field. If we succeed in coupling only one of the two zeromodes to a gauge 
field, we can hope to use this in order to construct a chiral gauge theory on the lattice. 
This appears to be impossible if one also insists that gauge invariance is maintained. 
However, if we do not worry about gauge invariance for the moment, we can couple 
the right handed mode to a gauge field, by replacing the free (2n-dimensional) Dirac 
operator and Wilson term by the gauge invariant ones, but only for a restricted number 
of s-slices around the right-handed domain wall, cf. [Q. In this way the gauge field 
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is confined within a "waveguide" around the domain wall, and interactions with the 
opposite chirality mode at the anti-domain wall are exponentially suppressed with Lg. 

We take the same gauge field on all s-slices inside the waveguide, which is natural if 
one thinks of s as a fiavor label, and define gauge transformations on the fermion field 
as _ _ 



WG = {s : So < s < s'o} (2.5) 

with in a gauge group G. The detailed choice of the boundaries Sq and s'q is not very 
important, provided they are sufficiently far from the domain wall that the zero mode is 
exponentially small at the waveguide boundary. For symmetry reasons we shall choose 
So = (-Ls + 2)/4 + l and Sq = (3Ls + 2)/4, such that the right handed mode at s = Ls/ 2 + 1 
is located at the center of the waveguide, see fig. 1. With this choice we have to take 
Ls — 2 a. multiple of four. 

Having made this division into a waveguide and its exterior, we note that the model 
has a global G x G symmetry: 

K - 9^1, % - seWG, 

■^i ^ /i^x, K ^ s^wG. ^ ' 

With our choice for the position of the waveguide boundary, there is a symmetry in- 
volving parity plus a refiection in the s-direction with respect to the plane s = Sq — \ = 
L,/4 + l, 

K - ld^pf^"-\ (2.7) 

with Px = (— Xi, ■ ■ ■ , — Xrf) the parity transform of x. 

It is clear that the hopping terms from sq ~ 1 to Sq and from s'q to s'q + 1 break the 
local gauge invariance of eq. ( |2.4|) . However, this can be repaired by putting in a scalar 
field V at the boundary of the waveguide, or alternatively by interpreting the gauge field 
Qx that appears in the action after performing a gauge transformation as a Stiickelberg 
field. This leads to the gauge invariant action 

seWG s^WG 

Y I^'Pl"^'^^ +^"^^Pr'^"]+J2^''^' (2-8) 

s^so-1,Sq * 
- 1/(F°"VPl^'° + FV^Pr^'"-^) - y(FV"fPL^'o+^ + F«^VPr^'o), 

where we have supplied the Yukawa term with a coupling constant y. Note that we take 
the same scalar field at both waveguide boundaries. Since we have chosen r = 1 we have 
written projectors in the hopping terms in s, Pr(l) = ^{^ + ■ PiP) and W{U) 



5 



are the usual gauge covariant Dirac operator and Wilson term, whose explicit form is 
not important here, since we shall only work with f/ = 1 in this paper. The field G G 
is the scalar field, which can be thought of as a (radially frozen) Higgs field, and which 
transforms as 

K ^ hV^gl (2.9) 
The transformation given in eq. ( p.7[ ) remains a symmetry if V transforms as 

V. - Vl,. (2.10) 

Since gauge invariance is broken in the model without scalar field, we add a mass 
term for the gauge boson, which on the lattice takes the form ^X)^ tr {U^ + f/^), with 
K the mass parameter in lattice units. It takes the form of a hopping term for V when 
this field is used to restore the gauge invariance of this mass term, 

Sv = -/tE tr (Kt/M-V;^^ + h.c). (2.11) 

X,fl 

Gauge invariance and the necessity to couple only the zeromode on one of the domain 
walls to the gauge field has led to an action which contains an additional scalar field. One 
might wonder whether there is a better way to introduce a gauge field which couples to 
only one of the domain wall zeromodes, but avoids the extra scalar field. Unfortunately, 
this appears to be difficult, if not impossible in a model which contains both domain 
walls, as we shall argue in sect. 6. For a proposal in a different direction, in which the 
anti-domain wall is avoided by keeping Lg strictly infinite, see ref. [ED . 



To get an idea about the physics of the model ( |2.8| ), we can start with ?/ = 0, in 
which case the scalar field is decoupled. However, now the gauged and ungauged parts 
of the action have decoupled completely as well, which implies that the two zeromodes 
on the domain walls are no longer balanced by each other. Therefore new zeromodes 
with opposite chirality must emerge which will be bound to the waveguide boundary. 
As an illustration we have plotted in fig. 1 the four zeromodes computed for the smallest 
plane wave momentum on a lattice with d = 2, Lg = 50 and U = 1. At y = (figure 
la) one recognizes the two expected massless modes at the domain walls, but also two 
modes at the waveguide boundary. These modes are masslessQ, because there can be no 
overlap between the left and right handed components across the waveguide boundary. 
For nonzero y the two components can overlap and they form a Dirac state with mass 
approximately equal to y, see fig. lb, where we took V = 1. One clearly sees how in this 
case the wave functions which are peaked at the waveguide boundary extend across this 
boundary. Note that the modes shown in fig. 1 are symmetric around the waveguide 
boundary, in accordance with the symmetry given in eq. ( p.7|) . The extra mirror modes 
at the waveguide boundary will be further discussed in the next section. 



^Of course, these modes are not exactly massless, because of the exponentially suppressed mixing 
between the domain wall and boundary modes. 
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2.3 Conjectured phase diagram 

To arrive at a chiral model, both this additional fermion at the waveguide boundary and 
the scalar field have to be decoupled. We first make the simplification of neglecting the 
gauge field dynamics by replacing [/ — > 1. This is reasonable, because we are interested 
in the scaling region at small gauge couphng. There we can write U/^x — ^if^^i^jxaj+zi 
with C/^ the gauge field in the smooth Landau gauge. The Q can be absorbed by a 
gauge transformation on i/j, i/j and V. Since is now smooth and close to one, we can 
treat this field in perturbation theory and put U = 1 in our numerical computations. 
Note that at large Yukawa couphngs the scalar field dynamics cannot be computed in 
perturbation theory. 

The light fermion modes are well localized, which implies that the mode at the do- 
main wall has only an exponentially small overlap with the scalar field at the waveguide 
boundary. With r = 1 the domain wall zeromode has a mag nitude oc e"™«^»/^ at the 
waveguide boundary, and for sufficiently large Lg the effective Yukawa coupling to the 
scalar field is exponentially suppressed, ye~'^°^''^'^. The new fermion mode, on the other 
hand, is localized at the s-slice that carries the scalar field and is coupled to it with 
strength y. Therefore we can use the freedom of adjusting the coupling constants y and 
K. to try to decouple the unwanted fields at the waveguide boundary, while keeping Lg 
sufficiently large as to ensure that the physics of the zeromode at the domain wall will 
remain unaffected. 

In the broken phase (or ferromagnetic (FM) phase), where the scalar field expectation 
value V = (Vx) is nonzero, we expect that the boundary fermion for small y gets a mass 
oc yv, but also the gauge boson acquires a mass oc v. Therefore we cannot decouple the 
fermion while keeping the gauge boson light. If we allow the gauge boson to be massive, 
it follows from the triviality of the Yukawa coupling in this region of the phase diagram, 
that the fermion mass will be of comparable magnitude. The remaining option is to 
choose K, in the symmetric phase (or paramagnetic (PM) phase). Here we expect from 
experience with the massive Yang-Mills model that the scalar field can be decoupled from 
the low energy physics of the fermion-gauge model, because deep inside the symmetric 
phase all scalar excitations will have masses of the order of the cutoff. However, since 
V — 0, one would also expect the boundary fermion to have mass zero and therefore not 
to decouple, which implies that the low energy model would be vectorlike. 

An interesting possibility is however that our model might exhibit the strong Yukawa 
coupling behavior found in other lattice Higgs- Yukawa models [25-27]. It was shown 
that for strong Yukawa couplings in the symmetric phase, such models exhibit another 
phase (denoted by PMS) in which the fermion and the scalar field form a massive bound 
state with mass of the order of the cutoff. We illustrate this desirable scenario with a 
possible phase diagram at /t = for our model, for the Yukawa couplings y and ye^"^oLs/4 
of the waveguide boundary and domain wall fermions respectively, shown in fig. 2. 

For L5 — > oo, the domain wall fermion has negligible Yukawa coupling, and we can 
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hope the phase diagram to be similar to that of the Yukawa model studied in ref. ^1 



for small y the system starts off in a weak symmetric phase (PMW); for increasing y 
the system comes into a broken phase, because the induced fermion interactions are of 
ferromagnetic nature. Then for still larger y the system enters a strong symmetric phase 
in which the boundary fermion becomes massive (denoted by PMSi). 

For finite values of Lg also the domain wall fermion gets strongly coupled for large 
y and this could take the system into a different symmetric phase (denoted by PMS2), 
in which both fermions are massive. Like between the PMW and PMSi phases, there 
may be a FM phase separating the PMSi from the PMS2, cf. fig. 2. We note in passing 
that the presence of two Yukawa couplings which are both proportional to y but differ 
by a large factor e"^"^"^'^ makes it very difficult to apply a strong coupling expansion in 
y, to investigate or establish the PMSi phase analytically. 

If the phase diagram of fig. 2 would be qualitatively correct for our model and a 
PMSi phase does exist, we could decouple the unwanted boundary fermion as well as the 
scalar field: we can choose the Yukawa coupling y sufficiently strong that the boundary 
fermion forms a bound state with the scalar field and acquires a mass of the order of the 
cutoff, whereas the domain wall fermion still is weakly coupled and remains massless. 
Sufficiently deep in this PMSi phase also the scalar field is very massive and should 
decouple. When we then turn on a smooth external gauge field inside the waveguide, 
the only light particle coupling to it is the right-handed fermion at the domain wall. For 
this scenario to work, we emphasize that the details of the conjectured phase diagram 
in fig. 2 are not important, but only that the PMSi phase exists. In the next section we 
shall investigate the scenario in more detail. 



3 Mirror fermion representation of the model 
3.1 Mode expansion 

The colloquial discussion in the previous section can be made more explicit by rewriting 
the action as follows. Relabel the right and left handed fermion fields, "^^l = Pr^l^^ 
as 

xi = XR = n-\ (3.1) 

and the same for r^l = ^Pl,r (note the reversal of L and R). The new label t runs 
from 1 to = Ls/2. In fig. 1 we have indicated this new labeling for the zeromode 
wave functions shown there. With our choice for sq, s'q and Lg we can define a domain 
wall mass for both fields and x, which is a step function in t satisfying, 

m* = m'o-^+* = m'»-*. (3.2) 



8 



With this relabehng the two domain wall zeromodes will reside in the Dirac fermion 
field ip, whereas the waveguide boundary zeromodes will reside in x- After substituting 
eq. into eq. ( |2.8|) with f/ = 1, the action turns into 



Lt 



t=i 

t=l t 
- yx\VPL + V^Pr)x' - yi^'-'iV^PL + V^P«)V^^'. 



(3.3) 



In this form, the action resembles that of an L^-fiavor mirror fermion model in the 
fashion of ref. with the fermion and x the mirror fermion field. In fact, for Lg = 2 



the hopping terms in t are absent, m* = and our model reduces to the mirror fermion 



model of ref. I^Sj with equal Yukawa couplings for the fermion and the mirror fermion, 
and a vanishing single-site mass term. For Lg > 2 our model has a more complicated 
mass matrix (i.e. nondiagonal couplings among the fiavors s or t) and if our model is 
going to be more successful in decoupling the mirror fermion than the traditional mirror 
fermion approach, it must come from this mass term. 

The mass matrix for the Lt fiavors in our model is not diagonal but this can be reme- 
died by more rewriting. First we expand the fermion fields in a plane wave basis, which 
diagonalizes the Dirac operator and Wilson term, = J2p^^^^'4'pj — I^p^'^^^'^p- 
Here X]p is a normalized sum over the momenta on the d dimensional lattice, X]p 1 = 1- 
Then we can write, 



"^^x = EE ij^lh^l + ^xlKxl + xliwp + m^)'4)l + ^l{wp + m*)x* 



t=i p 

Lt-l 



t=i t 

rLt 



E 



- y E {xI{Vp-,Pl + V,l,Pn)xl + %\ylpPL + ^p-.^kX* 

pq 



(3.4) 



with = 1(1 sin(p^), Wp the diagonal form of the Wilson term, Wp = I]p(l ~ cos(p^)) 
and Vp the Fourier transform of V^. For y = the action has the schematic form 



i>x 



x) 




(3.5) 



with M a. {p dependent) matrix in fiavor space, which can be read off from eq. (|3.4| ) 
This action can be diagonalized by making unitary transformations on ip and Xi 



LU-' 



ftX 



(3.6) 
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such that GfgMstFtg = fJ'fSfg. The matrices F and G are eigenfunctions of M^"M and 
MM''' respectively, labeled by the index /: 

{M^M),tFtf = l/i/pF,;, {MM^),tGtf = (3.7) 

For suitable choices of the phases of the eigenfunctions, the /i/'s are real. Substituting 
the mode expansion (|3.6|) into the action (|3.4|) with the momentum label restored, we 
arrive at 

/=i p 

-y E %G%{V,.,Pl + Vl,PR)G%il + uJl,F%^{Vl,PL + V,.,Pr)FI^uI . 
fg,PQ 

In this representation of the model, it is seen that all fermion modes u-^ and ^-^ 
interact with the scalar field, but that their effective Yukawa coupling is determined by 
the magnitude of their wave function at the waveguide boundaries t = 1 and t = Lf. For 
y = the model is seen to describe free, degenerate fermions and mirror fermions with 
momentum dependent mass /i^ (for /i^ ^ 0, the eigenstates are uj^ + and u)^ — ^0. 
Exactly one flavor, which we denote with / = 0, has /i^ = (up to terms exponentially 
suppressed in Lg) for \p\ < pc, where pc is the critical momentum defined in sect. 2.1. 
For r = 1 and tuq close to 1, the critical momentum is pc ~ V2- All other /i^ and also /i^ 
for p outside the critical momentum region, are 0(1) in lattice units. This is illustrated 
in fig. 3, where we show the lowest three masses as a function of the momentum (again 
we have chosen d = 2). 

This shows that for y = and momenta |p|^Pc, the model contains a massless 
fermion, uj^, as well as a massless mirror fermion, All other modes (/ ^ 0) as well 
as the species doublers have a mass of the order of the cutoff. The species doublers 
of the zeromode / = are massive because /ip is 0(1) for momenta with = ±7r. 
Furthermore, it is seen in fig. 3, that fi^ is almost exactly zero (it is exponentially small 
oc exp(— moi^s/4)) until it quickly rises to nonzero values for \p\ > pc- 

As was discussed already in sect. 2.1, fig. 1 shows the t-dependence of the zeromodes 
F*^ and G*^ of the fermion (indicated by in the figure) and mirror fermion (indicated 
by x) fo'^ the smallest momenta |j9| = tt / L ^ p^. It shows that the zeromode for the 
fermion is sharply peaked at t = (Lf + l)/2, i.e. at the domain wall and the zeromode for 
the mirror fermion is localized at the boundary, at t = L^. The non-zero modes, which 
are not shown in this figure, are not localized. 

3.2 Reduced model 

The action (|3.8| ) is an exact representation of the action for the domain wall fermions. 
The reason for writing it in this form is that it reveals, more clearly than the original 
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action, which fermion modes are important for the low energy physics. To shed hght on 
the model for ?/ 7^ 0, we shall exploit this separation of light and heavy modes in order to 
simplify the model by making a number of approximations, which we expect to hold for 
large Lg. First we shall neglect all non-zeromodes. This is a reasonable approximation, 
since these fermion modes have masses of the order of the cutoff, /i^ = 0(l/a), / 7^ 0. 
If they would couple strongly to the zeromodes, they could still be important, but from 
the Yukawa interaction in ( p.8|) one can see that such a coupling involves the overlap 
of a zeromode and a non-zeromode at t = 1 or t = L^. Since the non-zeromodes are 
not localized, the value of the wave functions \Gtf \ or \Ftf \ at any given t is of order 
1/\/LI and the contribution of Lg internal heavy flavor fermions is expected to be of 
order L^|Gi/|7/i^ = 0(1/ fif,). 

The remaining zeromodes have a momentum dependent Yukawa coupling. For the 
fermion 0;° this coupling is proportional to the squared absolute value of the wave 
function Fl^o which is exponentially small. Furthermore, the mixing with the mirror 
fermion is either exponentially small for momenta \p\ < pc, or the modes are very massive 
for large momenta \p\ > pc, and we discard such heavy modes in our approximation. 
Therefore in this approximation the model describes a free massless fermion with 
momentum cutoff at pc, and a mirror fermion with Yukawa coupling to the scalar field. 
This Yukawa coupling contains a momentum dependent factor GiqGIq, cf. eq. ( p.8|) . It 
turns out, however, that in the momentum range well below the cutoff Pc this factor is 
almost constant and close to one, and then quickly drops to a small value for \p\^Pc- 
This momentum dependence of jGiol is shown in fig. 4. This justifies the approximation 
that we also impose the momentum cutoff on the mirror fermion and neglect the wave 
function factor in the Yukawa coupling for \p\ < pc. 

All this leads to a simplified "reduced" model, described by the action 

s'^' = E H^P^'p + K^pO + y E lliv,l,PR + Vp-,PLK (3.9) 

|p|<Pc \p\,\q\<Pc 

Notice that this model differs from the mirror fermion model of ref. [2^ by the absence 
of a momentum dependent mixing term between fermions and mirror fermions and by 
the presence of the momentum cutoff \p\ < Pc- 

In this approximation, the model shows all the features discussed in the previous 
section. In particular we see that the zeromode a;° which in the full model is localized 
at the domain wall, is nicely decoupled from the mirror fermion which is localized at 
the boundary. The domain wall zeromode has an exponentially small interaction with 
the scalar field, which we neglected in the action ( p.9|) , but the mirror fermion couples 
to the scalar field with strength y. This mirror fermion will decouple if there exists a 
strong symmetric (PMSi) phase for large y in which the mirror fermion and scalar field 
form a massive bound state. To summarize, the action ( p.9|) should describe the physics 
of the full model for Lg 00, i.e. at the horizontal axis of the phasediagram in fig. 2. 
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The usual approach to show that such a phase exists is to write the action in terms 
of the fermion-scalar composite field, 



e; = T.iPR^,-<i + plv,-,)u (3.10) 

which is chosen such that the Yukawa term turns into a mass term for the fermion 
field ^' . This fermion does not transform under the gauge group G, hence we shall call 
it neutral. If the momentum cutoff were absent, we could invert this transformation, 
ip = J2q{PR^p-q + -PlV'JLp)^^, and a strong coupling approximation of the resulting 
action for C,' would predict a mass of the order of the cutoff for this neutral fermion. 
However, due to the momentum cutoff such an argument cannot be used for our model 
and in fact the transformation (|3.10|) leads to a nonlocal action for 



The momentum cutoff, which prevents a straightforward analytic demonstration 
that a strong coupling phase exists, makes this model markedly different from models 
which are known to have a strong symmetric phase. It is somewhat similar, however, 
to a fermion-Higgs model with hypercubical Yukawa coupling ^ . In these models 
the Yukawa interaction in momentum representation contains a momentum dependent 
factor, coming from the averaging of the scalar field over the hypercube, which suppresses 
the coupling strength for large values of the scalar field momentum. Such models are 
known not to have a strong symmetric phase. 

To summarize this section, we have shown that the domain wall fermion model can 
be rewritten as a mirror fermion model, with Lt = Ls/2 fiavors. In order to decouple 
the mirror partner of the domain wall zeromode, we must show the existence of a strong 
symmetric phase, where the mirror fermion forms a massive bound state with the scalar 
field. We have argued that the model for large Lg can be simplified to a reduced model 
with only a fermion and a mirror fermion. In this model we cannot show the existence 
of a strong phase using standard analytic techniques. The momentum dependence of 
the Yukawa interaction (which gives rise to the cutoff pc) is more similar to that of a 
fermion-Higgs models with a hypercubical Yukawa interaction than to models with a 
local Yukawa interaction. Models with a hypercubical Yukawa interaction that have 
previously been investigated, are known not to have a strong symmetric phase. 

Of course the similarity to hypercubically coupled fermion-Higgs models does not 
prove that a strong phase is absent in our model, and we shall search for it with numerical 
methods. The most direct way to show the existence of a strong coupling phase, is by 
measuring the mass of the boundary (mirror) fermion for strong Yukawa coupling. In 
the next section we shall study the fermion masses, both for the domain wall fermion 
and the boundary fermion, in the quenched approximation. We shall compare these 
results with the masses found from the reduced model. 
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4 Fermion spectrum: numerical results 



In order to substantiate the discussion in the previous sections, we shall compute the 
fermion spectrum of the full domain wall model in the quenched approximation. We 
expect from experience gained with other fermion-Higgs models, that the presence of a 
strong coupling symmetric phase if it exists, can already be shown within the quenched 
model. In the quenched approximation there are no real phase transitions separating 
the PMW, PMSi and PMS2 phases of fig. 2. One expects, however, that the FM phase 
separating these phases in the unquenched model, now turns into a cross-over region, 
which separates regions of the quenched phase diagram with different (weak and strong 
coupling) behavior. In the following we shall refer to these regions as weak and strong 
coupling phases, as in sect. (2.3). 

For weak Yukawa coupling the mirror fermion mass in the quenched approximation is 
expected to behave as ~ yv, where the scalar field expectation value v is zero in the 
symmetric phase and nonzero in the broken phase. A strong coupling symmetric phase 
would lead to mi;' = yc{K,), with c{k) a function of the scalar field hopping parameter 
K. Typically c(/t) > v{k,), it decreases with k, and in particular it is nonzero and 0(1) 
in the symmetric phase. For instance, in the model of ref. [Q, c{k,) ~ 1/z{k,) with 

oc tY{VxV^^^). For the domain wall fermion mass we expect nip ^ for all k, and 

For this numerical study we use the domain wall model in 2 + 1 dimensions with 
gauge group G=U{1), but we keep the gauge fields in the global symmetry limit U = 
1. The scalar field action (|2.11| ) then is that of an XY model, and in the quenched 
approximation, where the scalar field dynamics is determined solely by the action ( p.ll| ), 
there is a vortex phase and a spinwave phase. The Kosterlitz-Thouless phase transition 
is at K = fCc ~ 0.5 in our convention for the action. Of course spontaneous symmetry 
breaking does not really occur in this two dimensional model, but on a finite lattice the 
field expectation value v shows a behavior similar to that in a model with spontaneous 
symmetry breaking: it is nonzero and 0(1) for k, > and then quickly drops to a 
small (nonzero) value for k < k^. For increasing volumes v becomes closer to zero for 
K < Kc but also in the spinwave phase it decreases slowly, such that in the limit of 
infinite volume v = everywhere, as it should. We emphasize that in not too large 
volumes, in which there is a clear distinction between the value of v in the vortex and 
spinwave phases, we expect a similar relation between fermion mass and f as in a four 
dimensional model with spontaneous symmetry breaking. Hence we shall refer to the 
vortex phase as the symmetric phase and to the spinwave phase as the broken phase. 

To find the fermion masses, we have measured the propagator in momentum space, 

S'\p) = L-'J2e'^^--y\^l%), (4.1) 

with the two dimensional lattice volume. Optimally, one should measure the full 
matrix S''^*(p) in flavor space, for a number of small momenta p and from that compute 
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the massive and massless eigenstates. However, the number of flavors typically is large 
(we use for instance Ls = 26), and it is impractical to compute the propagator matrix 
for all flavors s. Since we are only interested in the masses of the light states and since 
we know that these states are localized either at the domain wall or at the waveguide 
boundary, it is sufficient to compute only S''^* for selected s- values s = 1, Sq — I, Sq, Ls/2 + 
1, Sq and Sq + 1. In fact, we know from the discussion in the previous section that the 
mirror fermion is localized near s = sq and we need only consider s = sq — 1 and sq if 
we are interested only in the mirror fermion mass. 

The parity symmetry of eqs. ( p.7| ) and ( |2.10| ) can be used, after averaging over the 
scalar field, to relate certain RR and LL components of the fermion propagator: 

SUp) = Stl^'-'^'^'^'-'iPp), (4.2) 

where Pp is the parity reflected two-momentum. We have used this relation to average 
over the appropriate RR and LL components, in order to increase statistics. For the 
RR or LL component of a free fermion propagator we expect 

S{p)rr{ll) = -iZpisin^pi) - {+)ism{p2))/C^sm^{p^) +m^), (4.3) 

where Zp is a wave function renormalization constant and mp is the mass. In flg. 5a 
we have plotted the inverse of the averaged RR and LL components, as a function of 
X^/i sin^(p^). We used a lattice of size L'^Lg = 12^26 with Yukawa coupling y = 0.5, 
at K = 0.5 near the phase transition and ttiq = 1.1. We used antiperiodic boundary 
conditions for the fermions in the t-direction. The data have been normalized such that 
the slope (determined from the flrst and second point) is one. The straight lines are 
fits to the data and the good quality of these fits shows that the fermions are (nearly) 
free. 

To see if the reduced model resembles the full model also in a quantitative way, we 
have computed the inverse mirror fermion propagator in this model. This result is shown 
in fig. 5b, again normalized to slope one. The normalized mirror fermion propagator is 
in good agreement with the one computed in the full model. 

As anticipated the RR component of the domain wall fermion at s = Ls/2 + 1 (and 
the LL component at s = 1) has zero mass. The mirror fermion modes at s = Sq — 1 (the 
RR component) and at sq (the LL component) have a small mass which is consistent 
with mp = yv. All other components are seen to have a mass of order one in lattice 
units. In the same fashion we have computed the fermion masses at other values of k 
and y. In all cases we found that the domain wall fermion remains massless. 

Of course the most interesting results are those for the mirror fermion mass at small 
K in the symmetric phase and at large values of y. Unfortunately, the data here are 
subject to large statistical fiuctuations. Even after averaging over 3000 scalar field 
configurations at k = 0.1 and y = 10 we found that the errorbars on the propagator are 
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comparable with the signal. The reason is that the propagator itself is very small, which 
prohibits a reliable analysis of the fermion propagator for such values of the couplings. 

Instead of using k, ^ such that f ^ 0, we can also look for strong coupling behavior 
at larger n. In the broken phase the weak and strong regions are less pronounced, 
but the presence of a nearby PMS phase should still show up in a deviation from the 
relation uif ~ yv. Since v decreases for k \ Kc, a characteristic feature of weak 
coupling behavior is a fermion mass which decreases As mentioned above, 

strong coupling behavior would show up through an opposite trend of the fermion mass 
as a function of k, increasing towards the phase transition. In fig. 6 we show the n 
dependence of the waveguide fermion mass at fixed y = 2. One recognizes the typical 
weak coupling behavior of the mass. From experience with other models we expect that 
y = 2 is already a strong coupling. For larger y it is difficult to measure the fermion 
mass reliably, because it is comparable to the cutoff in the range In fig. 7 we 

have plotted the y dependence of the mass for fixed k, = 0.5. Strong coupling behavior 
should show up as a relative increase of the mass compared to the weak coupling trend. 
From fig. 7, however, we can at most infer a relative decrease of the mass for y > 1. For 
comparison, we have also plotted the line yv{K = 0.5) in this figure. 

The crosses in figs. 6 and 7 are the masses obtained from the reduced model. One 
sees the same qualitative behavior as in the full model, but the masses are systematically 
higher (except when mp^2, which is beyond the cutoff, where, in the full model, mixing 
with all the other heavy modes presumably becomes important). This difference may 
be due to the momentum dependence of the fermion wave function in the full model. 
For increasing momentum and masses closer to the cutoff, we expect the wave function 
to spread out and the overlap at the waveguide boundary to decrease. This implies 
that the residue Zp of the fermion propagator ([4.3|) is not constant but decreases with 
increasing momentum. This also leads to an underestimate of the fermion mass in the 
full model. Keeping these systematic effects in mind, we consider the results of the 
reduced model in reasonably good agreement with the full model. 

Even though at this stage we do not yet find a conclusive answer for k deep in the 
symmetric phase and large Yukawa coupling, the results shown in figs. 5, 6 and 7, are 
consistent with the weak coupling mass relation mp = yv. Also the awkward behavior 
of the model at small k and large y is not what we expect from a model in the strong 
coupling phase. Only in the transition region between the two regimes we expect large 
statistical fiuctuations, but after the bound state has formed, the model should describe 
weakly coupled massive Dirac fermions, whose mass should be easy to measure. 

5 Search for a strong coupling symmetric phase 
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5.1 Eigenvalue spectra 

The resuhs for the mirror fermion mass described above are very suggestive but did not 
give a conclusive answer to the question whether a strong couphng phase exists in our 
model. Therefore we will attempt to approach this problem from a different angle in 
this section. The idea here will be that the presence of a strong coupling phase shows 
up in the distribution of the eigenvalues of the fermion matrix ||30|| . 

Of course, we would like to look at the eigenvalues of the domain wall fermion matrix 
directly. This ^LP'Lg x 2L^Ls (the factor 2 comes from the Dirac index) matrix M is 
obtained by writing the action ( |2.^ ), with U = 1, in the form S = ipMip. However, 
it is unpractical to study M directly, because this (nonhermitian) matrix is too large 
to handle numerically on reasonably sized lattices, and it is not clear what to expect 
for the distribution of the eigenvalues for M in the representation following from eq. 
( 12.81) . Only in the representation that diagonalizes the mass matrix in flavor space we 
should expect similarities with the eigenvalue spectra of free fermions (for small y) with 
momentum dependent masses. 

This suggests that we use the reduced model, which is formulated in terms of these 
mass eigenstates, and which contains much less degrees of freedom. The reasonable 
agreement of the results for the fermion masses discussed in the previous section supports 
this strategy. In the reduced model we can compute the distribution of the eigenvalues 
of the fermion matrix at small, intermediate and large Yukawa coupling. Then we can 
compare these eigenvalue spectra with those obtained in models for which a strong 
coupling phase is known to exist or to be absent. As such reference models we use a 
model with naive fermions with local Yukawa coupling, which has a strong coupling 
phase, and the same model with hypercubical Yukawa coupling, which has no strong 
coupling phase. The actions for these models, which we shall refer to as the Yic and 
models are. 

Sic = Y.^JxytPy + yY.^AV,PR+V:PL)^,, (5.1) 

xy X 

She = Y.'^Jxyi^y + yY.\T.'^AVx-bPR + V:_,PL)iJx. (5.2) 
xy X b 

The sum over h in the hypercubical Yukawa interaction runs over the four corners of the 
elementary plaquette, 6^ = 0, 1. After Fourier transforming the Yhe model, we find 

She = J2 '^p^^p'^P + yJ2 hp-gt/jpiVp-qPn + V*_pPL)iJq (5.3) 
p PI 

which contains a factor hp_q = H/^ e*^^''"'^''^/^ cos((p^ — g^)/2), which goes to zero for 
large momenta \Pf_i — — > vr of the V field. 

The results of this comparison are presented in fig. 8. Fig. 8a contains the spectra 
of our reduced domain wall model at k = 0.1 and y = 0.2, 1.0 and 4.0. We have plotted 
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the eigenvalues obtained from 5 quenched scalar field configurations, with lattice size 
L = 12. The figure shows that the scattering of the eigenvalues caused by the strongly 
fiuctuating scalar field increases with increasing y, as expected. However, there is no 
sign of a qualitative change for larger y. Also for Yukawa couplings ?/ > 4 we found that 
the spectra do not change qualitatively, they just scale proportionally to y. 

This can be contrasted to the y dependence of the eigenvalue spectra in the Yic model 
shown in fig. 8b. Here we see an increase of the fiuctuations for y growing from to 
1, then the eigenvalues A start to rearrange themselves along the boundary of a crude 
circle, which cuts the real axis at approximately ±y, such that the region around the 
origin becomes depleted of eigenvalues. This signals the existence of a strong coupling 
phase for y^l, cf. ref. [pO| . 

The spectra of our reduced model do not show such a qualitative change for large y 
and are much similar to those shown in fig. 8c, which were obtained from the Y/jc model, 
and which does not have a strong coupling phase. Since we expect the reduced model 
to be qualitatively similar to the full model for large L^, this result casts serious doubt 
on the existence of a strong coupling phase in our domain wall fermion model. 

The properties of the eigenvalue distribution of the fermion matrix are also refiected 
in the behavior of the conjugate gradient (CG) inversion. For small Yukawa coupling 
and using anti-periodic boundary conditions to regulate the zeromode for the fermions, 
we expect a rapid convergence of the CG inversion on our relatively small lattice. Then 
for increasing y the inversion rate should deteriorate, i.e. the number of CG iterations to 
reach the solution to a given precision will increase. If there is a strong coupling phase, 
the number of iterations reaches a maximum at the cross-over to the strong phase 
and then decreases again, because in the strong coupling phase the composite fermions 
are again weakly coupled and massive. In ref. it was found that the number of 
CG iterations provided an accurate indicator for the location of the cross-over and the 
existence of the strong coupling phase. 

In fig. 9 we have plotted the number of CG iterations as a function of ?/ at k = 0.1 
[L = 12, Ls = 26). One recognizes the expected rise of the number of iterations when 
y increases from to ~ 1.5. But, unlike what one expects for a model with a strong 
coupling phase, there is no decrease for large y. Also after y = 2 the number of iterations 
keeps rising, albeit at a slower rate and with larger fiuctuations than at small y. For 
comparison we have also plotted the y dependence of the number of CG iterations 
obtained in the Yic and Yhc model. In the Yic model, which has a strong coupling phase, 
the number of CG iterations clearly shows a peak at y ~ 1; in the Y^c model, which has 
no strong coupling region, we see a behavior similar to that of our domain wall model. 

We do not have an analytic method to establish the existence or nonexistence of 
a strong coupling phase in our model, but by comparing the Yic and Yhc models, one 
might conjecture that a strong phase can only exist if the fermion and scalar modes 
are coupled strongly over the full momentum range, including the high momenta modes 
with near ±7r. This is the case in the Y/c model, but both in the Yhc model and in 
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our domain wall model the Yukawa coupling is suppressed for large momenta, though 
the details of the momentum dependence are different in the two models. 

5.2 Dynamical fermions 

We have also attempted an unquenched simulation of the model. A direct simulation 
of the model with the action ( |2.8| ), S = ipMip, is not feasible, because DetM is not 
positive definite. Therefore we have simulated instead a model with an extra fermion 
field X added, with action S = x^^X^ such that the fermion determinant is given by 
Det(MtM). In this model we can use a hybrid Monte Carlo algorithm to include the 
fermions, because the fermion determinant is now manifestly positive definite. 

To look for a strong coupling symmetric phase, one should use a small value of k, 
such that the model is in the symmetric phase at y = 0. For y > we expect that the 
system switches to the broken phase. This disappearance of the symmetric phase for 
arbitrarily small but nonzero y is a special feature of two dimensional Yukawa models, 
cf. e.g. 1^ and references therein. Then for larger y we expect to find a symmetric 
phase if such a phase exists. Unfortunately, the system turns out to be extremely hard 
to simulate numerically for small k. For k, = and small y we could still measure 
nonzero field expectation values, but for increasing y we had to decrease the trajectory 
length progressively more, to unacceptably small values (e.g. at ?/ = 1 we had to use a 
step size dt = 0.01 with 10 steps per trajectory to maintain an acceptance rate larger 
than 75%, and at ?/ = 10 we had to use dt = 0.001). This results in huge autocorrelation 
and equilibrium times, which make a realistic simulation unfeasible for values of y in 
the interesting region. Presumably this is due to large fluctuations in the eigenvalues of 
M^M, as is suggested by the spectra of the reduced model shown in fig. 8. In fig. 9 we 
have also plotted the number of CG iterations required for inversions in the unquenched 
model at K = (full triangles). This shows the same steady increase with y as found in 
the quenched model. Also a tentative run at y = 10, showed none of the improvement 
we would expect after moving into a strong coupling phase. 

6 Other ways to couple the gauge fields 

The model we have studied in this paper resulted from an attempt to couple gauge fields 
to the chiral mode on only one of the domain walls, while preserving gauge invariance. 
This gave rise to an extra scalar field and additional mirror fermion modes at the 
waveguide boundary. These mirror fermion modes do not seem to decouple, in other 
words, there does not seem to exist a region in the phase diagram where the mirror 
fermions have masses of the order of the cutoff, while the zeromodes at the domain walls 
remain light. However, one might wonder whether the gauge field cannot be coupled 
to the fermions in a different way, which avoids these complications. Unfortunately, it 
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appears to be difficult to find a different approach without obvious flaws. Let us briefly 
discuss the original proposal |T2[, and its relation to the model studied in this paper. 



In the original proposal gauge fields were put on all links of the d + 1 dimensional 
lattice, i.e. the gauge field was taken to be a full, d + 1-dimensional gauge field. The 
action for these gauge fields was chosen as 

s{u) = E ( /5 E Re tr {u;,u:^^^u;l^M) + z^^+i E tr (f/;.K%^t/;^^V;t) 

(6.1) 

The gauge field in the d + 1 direction is denoted by V. By choosing the coupling P^+i 
for the extra field V sufficiently different from the plaquette coupling (5 for the gauge 
fields U^, it was hoped that at the domain wall the gauge field dynamics would still be 
d-dimensional at scales much below the cutoffi 

Eq. ( |6.1|) can be viewed as the action for a number of (i-dimensional gauge fields 
(labeled by s), coupled to equally many unitary scalar fields V^. For = 1 we have 
just Lg independent nonlinear sigma models in d dimensions, each with a critical point 
at (3d+i = Pc- For each s, the global symmetry group is G x G, with transforming as 

V' -^g'V^{g'+y, (6.2) 

with the in G. The full symmetry group Q = G^" is gauged by the ci- dimensional 
gauge fields The hopping terms in the s-direction in the fermionic part of the action 
look like the Yukawa terms in our model, eq. (|2.8|) (with y = 1): 

S^y = - E(^V^Pi^^+^ + "^'^W'^Pr^'). (6.3) 

s 

For U'^ = 1 and Pd+i < Pc, the symmetry Q is unbroken, and (V^"*) = for all s. It is 
then easy to see that in a mean field approximation, where is replaced by f = 0, the 
fermion action is that of 2 massless Wilson fermions and — 2 Wilson fermions with 
mass ~ mo, which are decoupled from each other, and vectorlike in d dimensions. This 



has been investigated in more detail in ref. [22 



If we now take Pd+i > (3c, the group G breaks down to its diagonal subgroup, 
G^" —>■ G, and only one gauge field remains massless. The other gauge fields would get 
a mass oc f = {V), and could be made very massive by choosing v at the cutoff. The 
fermion hopping terms in the s direction would survive (in mean field), and we would 
find the usual zeromodes at both domain walls. However, the massless gauge field is 
independent of s and couples equally to the modes at the domain and anti-domain walls, 
again rendering the model vectorlike. 

In a sense then, the model which we studied in this paper, is an improvement on this 
situation. Formally, our model corresponds to choosing P^+i = oo for all but two s-slices 
(where we set Pd+i = f^) , forcing the d + 1-dimensional gauge field to be d-dimensional. 
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There is no interaction with the antidomain wall if we choose the plaquette coupling 
P = CO outside the waveguide. 

In our model we have chosen an s-dependence of (3d+i and f3, (and of the fermion 
hopping parameter in the s-direction by the introduction of a Yukawa coupling y ^ 1 
at the waveguide boundary), which we think had the best chance of producing a chiral 
model. Of course, a more general s-dependence is possible, but we do not believe that 
this will improve the situation as described in this paper. 

7 Summary and conclusion 

In this paper we have considered a gauge theory with domain wall fermions in a fi- 
nite volume with a right-handed zeromode living at the domain wall and a left-handed 
zeromode at the anti-domain wall. The right-handed mode at the domain wall is cou- 
pled to a four dimensional gauge field which is confined to a waveguide around this 
domain wall. The left-handed mode at the anti-domain wall remains uncoupled p^ . 
The fermion hopping terms across the waveguide boundaries break gauge invariance, 
which is restored by promoting these hopping terms to Yukawa couplings in a way sim- 
ilar to the way fermion mass terms are made gauge invariant in the Standard Model. 
This leads to the introduction of a scalar field which lives only at the boundaries of the 
waveguide. There are two parameters in this model associated with this scalar field, a 
Yukawa coupling y and a hopping parameter k, (or equivalently a mass) for the scalar 
field. In our numerical work we have studied the scalar-fermion dynamics in the model 
with U(l) gauge symmetry in 2 -|- 1 dimensions. The gauge fields, which can be treated 
perturbatively, are switched off and we mainly used the quenched approximation. 

For vanishing Yukawa coupling the regions inside and outside the waveguide decouple 
from each other and from the scalar field. Therefore one would expect that new chiral 
zeromodes show up at the waveguide boundaries. This is indeed what happens: there is 
a left-handed mirror mode just on the inside of one of the waveguide boundaries, and a 
right-handed mirror mode on the outside (cf. fig. 1). The inside mirror fermion couples 
to the gauge field in the waveguide, resulting in a vectorlike theory. To show explicitly 
that our domain wall fermion model can be interpreted as a mirror fermion model, one 
can view the extra dimension as a flavor space. The hopping and single site terms in the 
extra dimension then generate a mass matrix, which is not diagonal in flavor space. By 
diagonalizing this mass matrix at y = 0, one recovers the massless domain wall modes 
as well as the massless mirror partners at the waveguide boundary. All other modes 
have masses of the order of the cutoff. For y ^ all modes have Yukawa interactions 
with the scalar field, proportional to y and to the magnitude of the wave function of 
the particular mode at the waveguide boundary. Since the wave function of the domain 
wall mode is exponentially small at the waveguide boundary, its Yukawa interaction is 
very weak, even at large values of y; the mirror mode, however, interacts strongly with 
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the scalar field. 

The crucial question is then whether the mirror fermion can be decoupled. A fa- 
vorable possibility would be that for large Yukawa coupling, the mirror fermion at the 
waveguide boundary forms a bound state with the scalar field, with a Dirac mass of the 
order of the cutoff, while the gauge symmetry remains unbroken. Such a strong cou- 
pling behavior is known to exist in many fermion- scalar models. In particular a strong 
symmetric or paramagnetic (PMS) phase has been established in these models. The 
key point in our model is that only the mirror fermion should become heavy, while the 
modes at the domain walls should remain massless. This would be conceivable, because 
the mirror mode couples much more strongly to the scalar field than the domain wall 
mode. 

One might ask whether such a scenario is excluded by the simple consideration that 
the massless mirror mode is required to cancel the anomaly generated by the domain 
wall mode. We think that this is not the case. If one turns on a smooth external 
gauge field, a Goldstone-Wilczek current will carry charge away from the domain wall. 
However, since no gauge field is present outside the waveguide, this current vanishes in 
that region, and the charge will have to be deposited somehow at the waveguide wall. 
Of course, if massless mirror fermions are present, they will do the job, much as the 
antidomain wall zeromodes did in the case without a waveguide, but rather with an 
external gauge field present throughout space-time [O, ITR ITH]. However, an alternative 



possibility is that a Wess-Zumino current carries the charge at the waveguide without 
any massless fermion modes being present. In the V = 1 gauge (where V is the scalar 
field), a charge density of the form jo oc eoijdiAj can be nonzero due to the discontinuity 
of the gauge potential A at the waveguide boundary. 

In our numerical work we have used several approaches to search for the existence 
of a PMS phase, but never found an indication that it exists. Our best evidence that it 
is absent, comes from computations of the mirror fermion mass at values of k near the 
phase transition to the symmetric phase, where the dependence of the mirror fermion 
mass on k, and y is consistent with the weak coupling mass relation mp = yv even for 
large values of y {v is the scalar vacuum expectation value). No sign of the behavior 
typical of a strong Yukawa coupling region was found (cf. figs. 6 and 7). For all values 
of the Yukawa coupling that we have considered, the domain wall zeromodes remain 
massless and unaffected by the Yukawa interactions. It would of course be nice to 
directly measure fermion masses deep in the symmetric phase, for small k and large y. 
However, in this parameter region the signal for the propagator disappears in the noise 
and we have not been able to obtain data with small enough errors to draw any definite 
conclusion about the mirror fermion mass. 

Since direct computations of the mirror fermion mass run into numerical difficulties, 
we have also studied a "reduced" model, which contains only the mirror fermion inter- 
acting with the scalar field. It is obtained from the full model for large Lg (the extent 
in the extra dimension) by discarding all fermion modes with masses of the order of 
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the cutoff. Since the Yukawa couphng of the domain wall zeromode is exponentially 
suppressed in L^, this mode can also be discarded. The momentum dependence of the 
mirror fermion wave function is such that the effective Yukawa coupling for this mode is 
suppressed for momenta larger than the critical momentum pc and hence we also discard 
these large momentum modes. We computed the eigenvalue spectrum of the fermion 
matrix for this reduced model as a function of y, and compared this with typical eigen- 
value spectra for simple Yukawa models. The eigenvalue spectra of such models are very 
different for large y, depending on whether a PMS phase does or does not exist pO| . 
The eigenvalue spectrum of the reduced model shows no sign of a PMS phase (cf. fig. 
8). 

Of course, one would like to study the eigenvalue spectrum of the full model directly. 
This was not possible due to the prohibitive amount of computer resources that would be 
needed. We beheve however, that the "reduced" model captures the essential features of 
the full model, one of which is the existence of an effective momentum cutoff at Pc- This 
belief is supported by a reasonable agreement between light fermion masses computed 
in the full and reduced models. 

It appears that the existence of an effective momentum cutoff pc in the theory, is 
the underlying reason for the failure to find a PMS phase. In the domain wall approach 
the fermion doublers are decoupled by making them heavy which implies that these 
modes are not bound to the domain wall or waveguide boundary, as is the case with the 
light modes. This implies that for large momenta near the doubler momenta = vr, 
the wave function of the boundary mode will be spread out in the extra dimension 
and it will be small at the location of the scalar field. Therefore the effective Yukawa 
coupling for these large momenta modes is necessarily small. However, this suppression 
of the Yukawa coupling for large momenta then prevents the formation of fermion- scalar 
field bound states necessary to have a strong coupling phase. As a result, the mirror 
fermion at the waveguide boundary stays light, and renders the theory vectorlike in the 
scaling region. If this picture is right, it points at a fundamental problem for domain 
wall fermions with a waveguide, not just for the two dimensional quenched U(l) model 
investigated here. 

The results discussed above were obtained in the quenched approximation. The 
unquenched model, assuming the decoupling of the boundary fermion would have been 
successful, would describe a single right handed fermion interacting with a U(l) gauge 
field. This model is anomalous and one might fear that our unfavorable results are a 
reflection thereof. This, however, is not likely, because we can also think of our model as 
the quenched approximation of a vectorlike model, obtained by adding an extra mirror 
fermion: writing the original action ( p.8|) as 5" = ipMil), we can add an extra fermion 
field X with action S = x^'^X- For the additional x fermion, the handedness of the 
zeromodes at the domain wall and waveguide boundary is reversed and the model is now 
anomaly free. In the quenched approximation, however, the extra fermion is irrelevant 
and this model reduces to the one studied here. 
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We have performed some unquenched simulations in the model with the extra fermion 
included, using a hybrid Monte Carlo algorithm. The results are inconclusive due to the 
very large autocorrelation and equilibration times, but do not contradict the conclusions 
described above. 

All our numerical computations have been carried out within a restricted range of 
Yukawa couplings, as typically the signal to noise ratio deteriorated prohibitively for 
large values of y in the symmetric phase. Therefore, it is not logically excluded that 
some PMS hke behavior might be found at values of y beyond y 10 or so. In particular, 
we have not tried to investigate the existence of a PMS2 phase as described in section 
2.3. This phase would not be interesting, however, for the construction of a chiral gauge 
theory, since also the domain wall zeromodes would be strongly coupled to the scalar 
field and neutral with respect to the gauge charge. 

To summarize, we believe that all the evidence presented in this paper — the close 
resemblance to a mirror fermion model with hypercubical Yukawa interaction, the k and 
y dependence of the mirror fermion mass where we could measure it, the distribution 
of the eigenvalues in the reduced model for large y and the behavior of the conjugate 
gradient algorithm — indicates that a PMS phase does not exist in our model and that 
the mirror fermions, which exist as a consequence of the introduction of a waveguide, 
cannot be made heavy. Therefore, a vectorlike gauge theory will result when gauge 
interactions are turned on. In view of the discussion in section 6, we expect that this 
negative result is quite general for domain wall fermion models in which the volume in 
the extra dimension is kept finite at any stage in the definition of the model. 
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Figure captions 

Fig. 1: Wave functions of the four lightest modes with momentum \p\ = tt/L, on a 
2+1 dimensional lattice with L = 18, Lg = 50 and mo = 1.1. The solid (dotted) 
lines represent right (left) handed components. The ip and x ^i-re the fermion and 
mirror fermion located at the domain wall (vertical bar) and waveguide boundary 
(dashed band) respectively. We indicate both the s and t labeling defined in eq. 
(13.11). Figure a is for y = 0, figure b for y=0.5. 

Fig. 2: Sketch of a phase diagram which would make the the domain wall fermion 
model successful. The y and ye'^"^"/^ indicate the effective Yukawa couplings 
for the domain wall and waveguide fermion respectively and k = 0. The various 
phases are explained in sect. 2.3. 

Fig. 3: Momentum dependence of the three lowest mass eigenvalues fi^ at y = 0, 
obtained on a 2 + 1 dimensional lattice with L = 16, = 50 and mg = 1.1. The 
domain wall fermion and waveguide fermion have degenerate masses. 

Fig. 4: Momentum dependence of the three wave functions corresponding to the 
eigenvalues fi^ shown in fig. 3, evaluated at the waveguide, t = 1, on a 2 + 1 
dimensional lattice with L = 16 and Lg = 50 and mo = 1.1. The symbols 
correspond to those of fig. 3. 

Fig. 5: Inverse propagator S^^{p) measured at s = 1, Sq — 1, Sq, Ls/2 + 1, s'q and Sq + I, 
for y = 0.5 and k = 0.5 on a 12^26 lattice with mo = 1.1. We show averaged 
components as explained in the text, errorbars are smaller that the symbols. For 
free naive fermions the fits with Ansatz S~^{p) = {rn?p + "^^^siv? p^) / Zp (solid 
lines) would be exact. Figure a is for the full model which has the domain wall 
zeromode (squares), the light waveguide mode (triangles) and many heavy modes 
(circles); figure b is for the reduced model, which only contains the waveguide 
mode. 

Fig. 6: The k dependence of the waveguide fermion mass at strong coupling y = 2, 
on a 12^26 lattice with mo = 1.1. The boxes (crosses) are for the full (reduced) 
model. 

Fig. 7: The y dependence of the waveguide fermion mass near the phase transition at 
K, = 0.5, on a 12^26 lattice with mo = 1.1. The boxes (crosses) are for the full 
(reduced) model. 

Fig. 8: Eigenvalue spectra for the reduced domain wall fermion model (figures a), the 
reference Yukawa model with local (figures b) and hypercubical (figures c) cou- 
pling. The left, middle and right figures are for y = 0.2, 1.0 and 4.0 respectively. 
The lattice size is = 12^ and k = 0.1. 
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Fig. 9: Number of conjugate gradient iterations to reduce the residual to < 10^^. The 
full symbols are for the quenched (circles) and unquenched (triangles) domain 
wall fermion model on a 12^26 lattice. The open symbols are for the Yukawa 
models defined in eq. ( p. 21 ) with local (triangles) and hypercubical (boxes) Yukawa 
interaction, on a 12^ lattice with k = 0.1. 
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